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Abstract: We calculate the pressure of massless (/>^-theory to order g^log{g) at weak 
coupling. The contributions to the pressure arise from the hard momentum scale of order 
T and the soft momentum scale of order gT. Effective field theory methods and dimensional 
reduction are used to separate the contributions from the two momentum scales: The hard 
contribution can be calculated as a power series in g'^ using naive perturbation theory with 
bare propagators. The soft contribution can be calculated using an effective theory in three 
dimensions, whose coefficients are power series in g^. This contribution is a power series in 
g starting at order 5^. The calculation of the hard part to order g^ involves a complicated 
four-loop sum-integral that was recently calculated by Gynther, Laine, Schroder, Torrero, 
and Vuorinen. The calculation of the soft part requires calculating the mass parameter in 
the effective theory to order g^ and the evaluation of five-loop vacuum diagrams in three 
dimensions. This gives the free energy correct up to order g"^ . The coefficients of the 
effective theory satisfy a set of renormalization group equations that can be used to sum 
up leading and subleading logarithms of T/ gT. We use the solutions to these equations to 
obtain a result for the free energy which is correct to order g^ \og[g). Finally, we investigate 
the convergence of the perturbative series. 



Keywords: 



Thermal Field Theory, NLO Computations, Renormalization Group , 



Contents 



|l]. Introduction 



|2|. Effective field theory 

.1| Coupling constant 

Coefficient of unit operator 
Mass parameter 



|3|. Soft contributions 
^ Results and discussion 
|5|. Summary 
\A[ Sum- integrals 



A.l One- loop sum-integrals 



A. 2 Two-loop sum-integrals 



A.3| Three-loop sum-integrals 



A.4 Four-loop sum-integrals 



Three-dimensional integrals 



B.l One-loop integrals 



B.2 Two- loop integrals 





n 

15 
19 
20 

m 



22 

n 
m 



B.3 Three- loop integrals 



B.4 Four-loop integrals 
|B.5| Five-loop integrals 



in. Explicit calculations 



25 



1. Introduction 

In recent years there has been significant progress in our understanding of thermal field 
theories in equilibrium ^, |3|, ^ . Part of the progress is based on the developement of the 
calculational technology necessary to perform loop calculations beyond the first correction. 
The motivation to carry out such difficult higher-order calculations of e.g. the pressure in 
thermal QCD is its relevance to heavy-ion collisions and the early universe. The pressure in 
nonabelian gauge theories has been calculated perturbatively through order g"^ in Ref. fs], 
to order in Refs. 0, §, and to order g^log{g) in Ref. @. There are three momentum 
scales that contribute to the pressure in thermal QCD - hard momenta of order T, soft 
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momenta of order gT, and supersoft momenta of order g^T. The next order - order g^ - is 
the first order at which all three momentum scales contribute to the pressure and it is also 
the order at which perturbation theory breaks down due to infrared divergences |10, 11 1. 
The pressure contains a nonperturbative contribution from the supersoft scale that can 
be estimated numerically [l^, |l^]. It also contains a presently unknown contribution 
from the hard scale. This contribution can be calculated by evaluating highly nontrivial 
four-loop vacuum diagrams with unresummed propagators. As a step in this direction, 
Gynther, Laine, Schroder, Torrero, and Vuorinen considered the simpler problem of (j)'^- 
theory at finite temperature and calculated the free energy to order g^ [|l^]. A difficult part 
of the calculation was to evaluate the four-loop triangle sum-integral, using the techniques 
developed by Arnold and Zhai in Refs. ^. 

In hot field theories at weak coupling, the momentum scales in the plasma are well 
separated and it is advantageous to use effective field theory methods to organize the 
calculations of the pressure into separate contributions from the hard, soft and supersoft 
scales. The basic idea is that the mass of the nonzero Matsubara modes are of order T 
and heavy. Since these modes are heavy, they decouple from the light modes, i.e. the 
static Matsubara modes. In particular, all fermionic modes decouple since their masses are 
always of order T. The contributions from the nonzero Matsubara modes to thermodynamic 
quantities can be calculated using bare propagators and are encoded in the parameters of 
the effective theory. Integrating out the hard scale T, i.e. integrating out the nonzero 
Matsubara frequencies, leaves us with an effective dimensionally reduced theory for the 
scales gT and g'^T |p. In the case of QCD, the effective theory is an SU{N) gauge theory 



coupled to an adjoint Higgs. The process is known as dimensional reduction [Qq, 17 



T|, H, H. The next step is to construct a second effective theory for the scale g^T by 
integrating out the scale gT from the problem Q. It amounts to integrating out the adjoint 
Higgs and this step can also be made in perturbation theory. This effective theory is a 
nonabelian gauge theory in three dimensions, which is confining with a nonperturbative 
mass gap of order g^T ||ll|. This theory must be treated nonperturbatively and gives the 
nonperturbative contribution to the pressure mentioned above. 

In the present paper we consider the thermodynamics of massless (/)^-theory and cal- 
culate the pressure through order g^ log(g') in a weak-coupling expansion using effective 
field theory. Calculations in scalar field theory are simplified by the fact that the supersoft 
scale g^T does not appear and so we only need to construct a single effective theory for the 
soft scale gT. This theory is infrared safe to all orders in perturbation theory due to the 
generation of a thermal mass of order gT. Compared to the (/^-calculations of Ref. |15|, 
the next order requires the matching of the mass parameter to three loops and the evalu- 
ation of some five-loop vacuum diagrams in the effective theory. The matching involves a 
nontrivial three-loop sum-integral that was calculated recently in Ref. [21|. 

The paper is organized as follows. In Sec. II, we briefiy discuss effective field theory 
and determine the coefficients of the dimensionally reduced theory. In Sec. Ill, we use the 
effective theory and calculate the soft contributions to the pressure. In Sec. IV, we present 
and discuss our final results for the pressure. In Sec. V, we summarize. In Appendix A and 
B, we list the necessary sum-integrals and integrals. In Appendix C, we calculate explicitly 



- 2 - 



some of the new three-dimensional integrals that we need. 



2. Effective field theory 

In this section, we briefly discuss the three-dimensional effective field theory and the 
matching procedure used to determine its coefficients. For a detailed discussion, see e.g. 
Refs. [11,1^. 

The Euclidean Lagrangian density for a massless scalar field with a ^^-interaction is 



where g is the coupling constant and A£ includes counterterms. This term reads 



(2.1) 



(2.2) 



In the present case we need the counterterm Ag'^ to next-to-leading order i g^. It is given 

by 



Ag' 



3 / 9 17 ^ 2 
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<7^ 



(2.3) 



where a = g'^ / (Att)'^ . We denote by (j){x) the field in the effective theory. It can be 
approximately, i.e. up to field redefinitions, be identified with zero-frequency mode of the 
field $ in the original theory. The Lagrangian of the effective theory can be then be written 
as 



/:efr = ^(V0)2 + imV + f|<^^ + 



(2.4) 



where m is the mass of the theory and g^ is the quartic coupling. The dots indicate an 
infinite series of higher-order operators consistent with the symmetries, such as rotational 
invariance and the discrete symmetry cp ^ —cp. In Eq. ( |2.4| ), we have omitted a coefficient 
/ of the unit operator. Its interpretation is that it gives the contribution to the free energy 
from the hard scale T. 

For the calculation of the pressure to order g^log{g), we need to know / and the 
mass parameter to order g^ and the coupling constant g^ to order g'^, i.e. we consider 
(/)^-theory in three spatial dimensions ^. This theory is superrenormalizable and only the 



mass needs renormalization |22]. The parameters in the effective Lagrangian (2^) are 
determined by calculating static correlation functions in the two theories at long distances 
R, i.e. R ^ 1/7", and demanding that they be the same [19|. In the matching calculations. 



we are employing strict perturbation theory [19|. This amounts to doing perturbative 
calculations in power series in g^ in which we treat the mass parameter as a perturbation 



""^ Power counting tells one that the operator {4>\7(j>)^ contributes to the free energy first at order 
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in the effective theory. The Lagrangian is therefore spht into a free and an interacting part 
according to 



/'free 



-eff 



/■int 
''-efr 



2 ^ ^ 24^ 



1 



(2.5) 
(2.6) 



Strict perturbation theory gives rise to infrared divergences in the calculation that phys- 
ically are cut off by the generation of a thermal mass m. The same infrared divergences 
appear in the loops in the full theory and so they cancel in the matching calculations. 
The incorrect treatment of the infrared divergences and the physics on the scale gT is not 
problematic since this will be taken care of by calculations in the effective theory. The 
matching calculations treat the physics on the hard scale correctly and the physics on that 
scale is encoded in the parameters of the three-dimensional effective Lagrangian. 

However, the matching calculations of the parameters in Ces are complicated by ul- 
traviolet divergences. Those divergences that are associated with the full four-dimensional 
theory are removed by renormalization of the coupling constant g. The remaining di- 
vergences are cancelled by the extra counterterms that are determined by the ultraviolet 
divergences in the effective theory. These divergences are regulated by introducing a cutoff 
A. The cutoff A can be thought of as an arbitrary factorization scale that separates the 



scale T from the scale gT (or smaller) which can be treated in the effective theory |19|. 
The parameters in the effective theory therefore depend on the cutoff A in order to cancel 
the A-dependence of the loop integrals in the effective theory. 



2.1 Coupling constant 

To leading order in the coupling g'^, we can simply read off 
the coupling g^ from the Lagrangian of the full theory. Making 
the replacement <I> — > VTcp in the Lagrangian (2.1) and comparing 
Jq dr C with £efri we conclude that g^ = g'^T. The one-loop graph 
needed for the matching of the coupling g^ to next-to-leading order 
in g^ is shown in Fig. Since the loop correction vanishes in 
the effective theory due to the fact that we are using massless 
propagators, the matching equation reduces to 



g^T 



(2.7) 



XX 



Figure 1: One-loop 
Feynman graph that 
contributes to the cou- 
pling g| in the effective 
theory. 



where Ai^f^ is the order-^f"^ coupling constant counterterm in Eq. ( |2.3| ). After renormaliza- 
tion, we find 



glik) = g\f,)T 



(47r 



1 



(47r 



AttT 



^(log^^ + 2,,log^ + "" 



AttT 



AttT 
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(2.8) 
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where = g^{n) is the couphng constant at the scale ^ in the MS scheme and we have 
kept the order-e terms in for later use. We have used the renormalization group equation 
for the running coupling constant (7^, 

5 9 17 o , , 

ijL—a = Sa^ , (2.9) 

Ofi 3 

to change the scale from A to ^. The right-hand side of Eq. ( p.^ ) is independent of A. In 
fact, since the coupling g'^ does not require renormalization in three dimensions, it satisfies 
the renormalization group equation 

a|^53 = . (2.10) 

2.2 Coefficient of unit operator 

The partition function in the full theory is given by the path integral 

Z = J p$e-^''^^J~'^'^^ , (2.11) 

and the pressure is then given hy V = TlogZ /V , where V is the volume of the system. In 
terms of the effective theory, the partition function can be written as 

Z = e-f^ j P(/.e-/'^'^'^-ff . (2.12) 

The matching then yields 

log^ = -/F + log^eff , (2.13) 

where ^eff is the partitition function of the three-dimensional theory. Equivalently, we 
can write T = JT^ard + -^soft, where J^hard = fT and J^^oit = -TlogZ^s/V . Now since 
calculations in strict perturbation theory in the effective theory is carried out using bare 
propagators, there is no scale in the vacuum graphs. This implies that they vanish in 
dimensional regularization and that log^^eflf = 0. Eq. ( |2.13| ) then tells us that / is given by 
a strict loop expansion in four dimensions. 

The vacuum diagrams through four loops are shown in Figs. 

o 00 

Figure 2: One-loop vacuum diagram. Figure 3: Two-loop vacuum diagram. 

We can then write 

_ ^{h) (h) (h) [h] (h) (h) [h] (h) 

>hard - -^0 + -^1 + -^2a + -^2b + -^3a + -^3b + -^3c + -^3d + 

+ -^{^,g' + W) + 2 + ^ Ai5^ , (2.14) 

r \ 9 9 
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Figure 4: Three-loop vacuum diagrams. 



(6) (c) (d) 

Figure 5: Four-loop vacuum diagrams. 



where Aig'^ and A.2g'^ are the oider-g'^ and ordev-g^ coupling constant counterterms, re- 
spectively, given in Eq. ( |2.3D . The superscript h indicates that the expression gives the 
hard contribution to the free energy. The expressions for the diagrams are 
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{h) 



(h) 



•J c 



{h) 



2b 



J 1 



(h) 



3a 



J 1 



3b 



3c 



{h) 

3d 



1 2/^ 1 



off 



16^ V^P^j 2/^g4 ' 



1 



1 



2 



32^ VTpP2 



1 



1 

3 



48" Tpge VTpP2 



(2.15) 
(2.16) 
(2.17) 
(2.18) 

(2.19) 

(2.20) 
(2.21) 
(2.22) 



where the symbol J] is defined in Eq. ( A.l ) and the self-energy n(P) is defined in Eq. ( |A1^ ). 
The expressions for the sum-integrals are listed in Appendix A. After renormalization, the 
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final expression is |15| 

-^hard(^) = — 
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C(-3) 
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+ Oie: 



(2.23) 



where a = aifi), C{^^^^ = 48.7976, Cf^i^^gi^ = -25.7055, and C^i^^gie = 28.9250. We have 
used the renormalization group equation for to cliange the renormalization scale from A 
to fi. Note that the final results contains a pole in e. We cancel it by adding a counterterm 
T6f [p. The term 6f can be determined by calulating the ultraviolet divergences in the 
effective theory. The triangle diagram in three dimensions has a logarithmic ultraviolet 
divergence and the counterterm needed to cancel this divergence is given by 



Sf 



(2.24) 



1536(47r)4e ' 

If we express the counterterm in terms of the coupling g of the full theory, we must take 
into account that g^ multiplies a pole in e and it therefore picks up finite terms. These 
terms will be of order g^ and can be neglected in the present calculation ^. The coefficient 
/ satisfies the evolution equation 



A 9 , 



192(4^)4^^ • (^-^^^ 

This follows from the scale dependence of the triangle diagram in three dimensions and the 
fact that the A-dependence of / must cancel the scale dependence in the effective theory. 

2.3 Mass parameter 

The simplest way of determining the mass parameter is by matching the Debye or 



screening mass rriD in the full theory and in the effective theory [19|. The Debye mass 
is given by the pole of static propagator, i.e. by 



/ + n(po = 0,p) = 



P 



(2.26) 



^Note that minimal subtraction in the full theory and in the effective theory are not equivalent. The 
difference is the finite terms mentioned above Iql. 
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where n(pO)P) denotes the self-energy function. In the effective theory, the equation is 



/ + + Ucsip) = , = -ml) , 



(2.27) 



where Ilefj (p) is the self-energy in the effective theory. Since the self-energy in the full theory 
is expanded around p = 0, we should do to the same in the effective theory (see discussion 
below). The loop integrals are therefore evaluated at zero external momentum and since the 
matching is carried out using massless propagators there is no scale in the loop integrals. 
They therefore vanish in in dimensional regularization, i.e. nefj(O) = n^fj(O) = ... = 0. 



Using this fact and equating ( |2.26| ) and ( 2.27 ), we obtain 



m 



m 



2 3 
D 



(2.28) 



The diagrams that contribute to the self-energy n(P) through three loops are shown in 
Fig |. 



O 



(a) 




(b) 



e 



(c) 






(e) (f) (g) (h) 

Figure 6: Fcynman graphs that contribute to the self-energy through three loops. 

The self-energy n(P) is given by 



n(p) = iif\p) + n^"^(P) + n^"^(P) + 



+ 



9 



(P) 2 . 2^ U-i^HP) 



(2.29) 



^Note that we use the symbol "r^" to emphasize that the the mass parameter is equal to the Debye 
mass m|j only in strict perturbation theory. The interpretation is that m gives the contribution to the 
Debye mass from the hard scale T. 
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The expression for the various terms in the self-energy are given by 



mP)='^.%^. (2.30) 



1 



1 



r %RQ-'R-'{P + Q + Rf ' 
mP)-\H^.{4^^?)\ (2.33, 



■Jh),^. 1 «^ 1 fU 1 



2 



Since the leading-order solution to Eq. ( p.28 ) gives a value of p that is of the order gT, it is 



justified to expand the loop diagrams in a Taylor series around p = 0. We can then write 



Eq. (2.28) as 



ml = nr(0) + UriO) + nr(0)/ + n3(0) + ... , = -ml , (2.38) 



or ml = nf ^(0) -hn^''^(0) -hfif \o) - ni(0)n^(0) . We then need the two-loop self-energy 

4?(P)andn(J 



diagram n2b(.P) to order p^, while the three- loop self-energy diagrams flij^ (P) and I[^{P) 



can be evaluated at p = 0. This yields 



= i/rf -^[n((3)|2 . (2.41) 
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The sum-integrals needed are listed in Appendix A. After renormalization, we obtain 
1 
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where g = g{A) and C/ = —38.4672. The mass parameter through order g^ is known to 
order e but we only need it to order e'^. We notice that the mass parameter contains 
uncancelled poles in e. It is advantageous to write the mass term as a sum of a finite piece 
and a counterterm Am^, where 
1 



m^(A) 



9'{f^)T' X 
2 



24 

xi 1 + 



9 



(47r)2 
65^ 



4 log 



A 
4^ 



3 log 



(4vr) 



41og^ 



A 
4^ 



47rr 
2 ^ 



+ 2 - 7s + 2 



3^ 

-log' 
2 ^ 47rr 



C^(-i) 
C(-i) 

, 19 C'(-l) 



1 ^ , 

log h 



. , A 2851 
+47^; log + 



AttT 



4rr2 



g^T 



24(47r)2e 
5|(A) 



I c^(-i) 

C(-i) 

+21og2(27r) 
69' 



113 _ 
"72" 12 



864 
17 

7s 



95 
48 



7i 



119 
144 



7E 



1 



C(3) - 771 + 



^C'ball + 



+ 0(e) 



(2.43) 



(47r) 
65' 



log 



47rr 



+ 7E 



(4^; 



2 I log'73^ + 27i?log ^ 
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- 24(4vr)2e ' ^'-''^ 

where g = g{fi) and we have used Eq. ( |2.9D to change the renormalization scale from A to 
jj,. The term Am^ acts as a counterterm in the effective theory. In fact, the sunset diagram 
in three dimensions that contribute to the self-energy is logarithmically divergent, whose 
divergence exactly is given by the right-hand side of Eq. (2.44) p^]. The mass parameter 
m in three dimensions therefore satisfies the evolution equation 

1 9l 



6 (47r) 



(2.45) 
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In the remainder of the paper, we will use m instead of rfi for covenience. 
3. Soft contributions 

In this section, we calculate the soft contributions Vgoft to the pressure. This requires the 
calculations of vacuum diagrams in the effective theory (2.4) through five loops. In order 
to take into account the soft scale gT, we now include the mass term in the free part 
of the Lagrangian and only the quartic term in Eq. (|2.4| ) is treated as an interaction. The 
inclusion of the mass term in the propagators cuts off the infrared divergences that plagues 
naive perturbation theory in the full theory. 

The one-loop vacuum diagram is shown in Fig ^. Its contribution to the free energy is 
given by 



(3.1) 



where the superscript (s) indicates that the expression gives the soft contribution to the 
free energy. Using the expression in the Appendix B, we obtain 



' Utt 



(3.2) 



The two-loop vacuum diagram is shown in Fig ^. Its contribution to the free energy is 
given by 



1 



Using the expression in the Appendix B, we obtain 

gWT 



8{4tt) 



(3.3) 



(3.4) 



The three-loop vacuum diagrams are shown in Fig ^. The contribution to the free energy 
is given by 
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where Am is the mass counterterm (2.44) in the effective theory and 
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Using the expression in the Appendix B, we obtain 
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We note that all poles in e cancel as they must since there are no divergences from the 
hard part proportional to g^m. 

The four-loop vacuum diagrams are shown in Fig ^. The contribution to the free 
energy is given by 



dm? 



(3.9) 



where the expressions for the diagrams are 
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(3.12) 
(3.13) 



Using the expressions in the Appendix B, we obtain 



768(47r)4 
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1536(47r)4e 



(3.14) 



The pole in e in Eq. ( 3.14| ) arises from the triangle diagram in Eq. ( 3.13| ). This pole is 



cancelled by the counterterm in Eq. ( 2.24 ). 

The five-loop vacuum diagrams are shown in Fig ^. The contributions to the free 
energy are given by 
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Figure 7: Five-loop vacuum diagrams that contribute to the soft part of the free energy. 



where the expressions for the diagrams are 
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pqrst Q'^ + (P + q)^ + ''^ + "T-^ (P + r)^ + 
1111 

^ 

+ (p + s)2 + + (p + t)^ + m 

1 «„ /■ 1 /■ 1 1 1 



2 



2 



(3.23) 



144 Jp (p2 + m2)2 g2 _|_ j^2 ^2 _|_ j^2 _j_ q _|_ j.~)2 _|_ j^: 

. / , (3.24) 

32 Jpqrst q"^ + (P + + (P + J^)^ + + r)2 + + 

111 / X 

X . (3.25) 

(p + s)2 + (s + t)2 + + 777,2 ' 



Using the expressions in the Appendix B, we obtain 



is) 



9lT 



288m(47r)5 

A_ 1 
2m 4 



log2 -iL + i (1 - 8 log 2) log A _ 1^ _ l^r^ + ^,,2 2+ 

2m 



15 

64 



23 81 
+— log 2 + 6 log2 2 - 6 log 3 - — C(3) + 5Li2(i) + 9C4j 



(3.26) 



where C^j = 0.443166. We note that all poles in e cancel as they must since there are no 
divergences from the hard part proportional to gf/m. Adding Eqs. (|3.2|) , (|3.3|) , (^), 
(3.14), and ( |3.26| ) as well as the counterterm Eq. ( p. 24 ), we obtain the soft contribution to 
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the free energy through five loops 



T 
•I I 



0+1+2+3+4 



127r 8(47r)2 96(47r)2 



81ogA + 9_81og2 



+ 



768(47r)4 



A 



4(4 - vr^) log 4 + 16 log 2 - 42C(3) + 7r^(l + 2 log 2) 



2m 



288m(47r)5 

" o A \ , , , , A 15 3 o 9 o , 
log' — + - 1 - 81og 2 log — - — - -vr^ + -^2 log 2+ 
2m 4 2m 64 8 8 

23 81 
+— log 2 + 6 log2 2 - 6 log 3 - — C(3) + 5Li2(i) + 9C4j 



(3.27) 



Using the evolution equations for g| and m^, it easy to check that the free energy, Eq. ( 2.23| ) 



plus Eq. (3.27) is independent of the factorization scale A. 

By expanding the coupling ( |2.8| ) and the mass parameter vr? (2.43) to the appro- 
priate orders in the various terms in ( |3.27| ), we obtain the soft contribution through order 
. This yields 

•>' soft — 



90 



r 5^6 



3/2 



15 



-a 



15^/6 



-a 



5/2 



log 



[I 2 
InT ~ 3 



log a + C5 



15 , 
-16"^ 



-48 log ^ + 16^^ - 327i, - 84C(3) + 8 + 16 log ^+ 

2 A \ 

+161oga + 7r^ ( 2 + 121og2-41og- -41oga + 81og^ j 



+ 



225v6 7/9 

' 2 /^221 2 4 2 4C'(-1) 4^ \ fi 

4^ + + 3^^ - 3 3 - 3 CMT " 3 '°^"J ^""^4^^ 

[2 8 C'(-l) 52 8 2\ ^ 4,2 

+ \ — T IE H log - log a H log"^ a+ 

\15 45 C(-l) 45' 45 3y 45 



C7 



where the constants C5 and C7 are defined below. 



4. Results and discussion 



(3.28) 



The full pressure is given by minus the sum of Eq. ( 2.23 ) and Eq. ( |3.27| ). The strict 
weak-coupling result for the pressure through order g''^ is minus the sum of Eq. (|]2|) and 
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Eq. (|3^ . This yields 

V = Pidcal X 



5 5\/6 15 2 

4 3 4 



log 



AttT 



15^/6 



-a 



5/2 



45 



-a 



log' 



log 



47rr 
1 



■ log a + Cs 



/269 

47rr 3 lis" 



9 



1 



+ ^{4-TT^)loga + Ce 



+ 



225 v6 

H T^a'''' X 



log' 



+ 



/221 2 

' \- TilE 



4ttT \135 3 



77 log o - o 



2 4C'(-1) 4 



3 3 C(-l) 3 



- log a log 



47rr 



52 8 2\ 4,2 

— 7_E log - log a H log a + u? 

15 45 C(-l) 45 ' 45 3y 45 



where Pideai = 7r^T^/90 and where the constants C4 — C7 are 
C4 = 

5 



^.1 4 CV1) 2 C^(-3) 
45 3^-^ 3C(-1) 3C(-3) 



C5 



6 + 



2 2 2^-1) 
— loe ^ — - 

3 ^3 3 C(-l) 



+ 

, (4.1) 

(4.2) 
(4.3) 



1,^ 2m 2 103 1 , 
3(4-vr^)log- + — + 



25 n 57r^ 

H 7f H 

36 24 



g ^triangle 



vr . 4 511 

-j^2 triangle "I" 180^^ 



7r2 2, A75 

y7i. + vr'log2+f — 



:1E 



Ci-1) 2 /C^(-l) 

c(-i) ^3Vc(-i) 



+ 



5C"(-1) 2 C'(-3) 2267 
^9 C(-l) ~ 3^^ CFsy ~ '324 



C(3) 



(4.4) 



1457 1 
810 ^45 



2 ^ 749 56 
15^^ + 270^^ + 15^^ 



11 2^ 2 C'i-l) 

H ^ 

20 15 C(-l) 



16 , , , 16, 2. N 52 , 2 19C'(-1) 



+ - 



34 , 2 2 2 , n 4 , 2 o 28 

H log - + -vr^ log 2 H log^ 3 H 

15 3 5 45 15 



+ 

C(-l) 
45 '^C(-l) 

log' 2 - 



38 



+ 



C^(-i) 

45 VC(-l) 

8 1 o 8C'(-1), 2 97^,^, 16,. 97 o 16^ 

log 2 log 3 H ^ — flog C(3) H Li2(t) H Tf H (^41 , 

45 45 C(-1 3 54^^^ 9 ^^4^ ^ gg 5 , 



(4.5) 



where C^,-,^^^^^ 



-25.7055 and C^,^^^,^ 

C4 = 
C5 = 
Q = 
C7 = 



= 28.9250. The numerical values of C4 — C7 are 
1.09775 , (4.6) 



-0.0273205 , 
-6.5936 , 
-0.862 . 



(4.7) 
(4.8) 
(4.9) 
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Note that the A-dependence cancels in the result (^j). Using Eq. (|2.9| ) for the running of 
a, it is straightforward to check that the final result Eq. ( [4.1D is RG invariant up to higher- 
order corrections. The ordev-g^ result was obtained by Frenkel, Saa, and Taylor p3| , 
the order-fl'^ result by Parwani and Singh |2^], the order-(7^ log(5) result by Braaten and 
Nieto [19 1, and the oicdei-g^ result by Gynther et al |15|. The latter was later reproduced 
in Ref. [21 1 using screened perturbation theory [25, 27[ by taking the weak-coupling 
limit for the mass parameter, m = gT/ \/24. 

An expansion of the pressure in powers of g is given in Eq. (|4.1|). It is accurate up to 



corrections of order g^\og{g). A more accurate expression can be obtained by using the 
fact that our short-distance coefficients satisfy a set of evolution equations. The solutions 
to the evolution equations are 



5i(A) 

/(A) 



5i(2vrT) , 

, 7r2o§(27rr) A 



m^{K) = m^{2-KT) + 



192(47r)4 
6(47r)2 



log 



27rr 
A 



2tiT 



(4.10) 
(4.11) 

(4.12) 



If we substitute the short-distance coefficients (|4.10D and ( [4.12 ) into Eq. ( 3.27 ) and add 
the short-distance contribution ( 4.11] ), setting A = gT/^/2A everwhere, and expand the 
resulting expression in powers of we obtain the complete result for the pressure, which 
is correct up to order g^\og{g). The contributions to the free energy of order g^\og{g) 
come from ( 4.11| ) and from using ( 4.12| ) to expand the g^rn?T term in ( |3.27| ). This yields 



log{g) = (log 2 - 7E) (4 - vr^) log(<7) • (4.13) 

Moreover, using the solutions to the flow equations, we are summing up leading logarithms 
of the form g'^^'^^ \og^ {g) and e.g. subleading logarithms of the form (7^""'"^ log"(g'), where 
n = 2, 3, .... These terms are obtained by expanding out the m?T and g^mT terms in (|3.27| ), 
respectively. 

In Fig. |, we show the various loop orders of "Phard normalized to "Pideai to orders g , 
g'^, and g^, where "Phard is given by minus Eq. ( p. 231) ^. We have chosen fj, = 2ttT and 
A = 2ttT. We notice that the successive approximations are larger than the previous one. 
In Fig. ^, we show the weak-coupling expansion of Pgoft normalized to Pidoai to orders g^, 
ff^! 9^ 9^1 9^ : and g^ \og{g), where Psoft is given by minus the sum of Eqs. ( p. 281 ) and ( 4.13| ). 

In Fig. ^ we show the weak-coupling expansion of the pressure V given by ( [4.1| ) 
minus ([4.13| ) normalized to Pideai to orders 5^, g^ g^, g^ g^, (7^, and g^log{g) . The 
convergence properties of the successive approximations of the sum P = P'hard + Psoft 
clearly is better than the convergence properties of the successive approximations to Phard 
and P'soft separately. 

In Fig. 11, we plot the successive loop orders of minus Eq. ( p. 27 ) normalized to P'ideai- 
In the one- and two- loop approximations, we use the leading-order results for g^ and for m^. 



Note that we omit the pole in e in Eq. (2.23) in the plots of the hard part. 
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Figure 8: Hard contributions T'hard to the 
pressure V normalized to Pideai to order g'^, 
g^, and g^. 
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Figure 9: Soft contributions T'soft to the 
pressure V normalized to T'idcai to order g^, 
9^, ff^ 9^, 5^ and g^log(5). 
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Figure 10: Weak-coupling expansion of the 
pressure V normalized to Pidoai to order g^, 
ff^ 9^, £'^ ff^ and 9^\og{g). 
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Figure 11: Soft contributions Vsoh to 
the pressure V normalized to Pideai at one 
through five loops. 



At three and four loops, we use the leading-order result for and next-to-leading order 
result for m?. Finally, at five loops, we use the solutions to the evolution equations for g^, f, 
and m?. The renormalization scale is n = 2-kT and the factorization scale is A = gT /^/TA. 
These approximations represent a selective resummation of higher-order terms. Clearly, 
the convergence is better than the strict perturbative expansion. In particular, the three-, 
four-, and five-loop approximations are very close. 



In Fig. 12, we plot the successive loop orders of the the pressure which is given by the 
sum of minus Eq. ( 2.23| ) minus Eq. ( |3.27 ), and minus (|4.13|) , normalized to "Pideal > starting 
at two loops. We are using the same approximations for g"^ and rn? as in the previous plot. 
Again we notice that the convergence of P is better than Phard and Pgoft separately. In fact 
the convergence is very good as the 3-loop through 5-loop approximations are very close. 
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It is not surprising that a selective resummation improves the convergence of the series. 
This was also notice in screened perturbation theory [25, 27, 21 1. 



5. Summary 

In the present paper, we have calculated the 
pressure to order log (5) in massless (/>^-theory 
at weak coupling. The first step is the determi- 
nation of the coefficients in the dimensionally 
reduced effective field theory. This calculation 
encodes the physics of the hard scale T. The 
mass parameter was needed to order and in- 
volves a nontrivial three-loop sum-integral that 
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Figure 12: Successive approximations to 
the pressure V normalized to T'idoai at two 
through five loops., 



was recently calculated in Ref. |21|. The sec- 
ond step consists of using the effective theory 
to calculate the vacuum diagrams through five 
loops. All loop diagrams in the effective theory 
but one could be calculated analytically with 
dimensional regularization. This way of orga- 
nizing the calculations is more economical and 
efficient than resummed perturbation theory. 

The parameters of the effective theory, g^, f, and m^, satisfy a set of evolution equa- 
tions. The solutions of these equations show that the parameters depend explicitly on the 
renormalization scale. This dependence is necessary to cancel the dependence on the scale 
in the effective theory |19[ |. The fact the our final result for the pressure is independent of 
the renormalization scale is a nontrivial check of the calculations. Furthermore, by choos- 
ing A = gT /\/24 and using the solutions to the evolution equations, we were able to sum 
up leading logarithms of the form (^^n+s iog"'(g^) and e.g. subleading logarithms of the form 
g2n+5 iog"'(5), where n = 2, 3, .... as well as obtaining the coefficient of the g^ log{g) term. 

As pointed out in Ref. [15|, it would be advantageous to develop the machinery of 
calculating complicated multiloop sum-integrals in an automated fashion as has been done 
for Feynman diagrams at zero temperature. Perhaps such techniques could provide an- 
alytical expressions for the constants that today are known only numerically. This is 
necessary if one wants to tackle the formidable problem of calculating the hard part of the 
^^-contribution to the free energy of QCD. 
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A. Sum-integrals 



In the imaginary-time formalism for thermal field theory, the 4- momentum P = (Pq, p) is 
Euclidean with P^ = Pq +p2. The Euclidean energy po has discrete values: Pq = 2miT for 
bosons, where n is an integer. Loop diagrams involve sums over Pq and integrals over p. 
With dimensional regularization, the integral is generalized to d = 3— 2e spatial dimensions. 
We define the dimensionally regularized sum-integral by 



where 3 — 2e is the dimension of space and /i is an arbitrary momentum scale. The factor 
{e^ [A-kY is introduced so that, after minimal subtraction of the poles in e due to ultraviolet 
divergences, jjl coincides with the renormalization scale of the MS renormalization scheme. 



A.l One- loop sum- integrals 

The massless one-loop sum-integral is given by 




(A.l) 





C(2n - 3 + 2e) r(n - I + e) 

8^^^ r(i)r(n) 



(27rT) 



4-2n-2e 



(A.2) 



where C,{x) is Riemann's zeta function. Specifically, we need the sum-integrals 




[1 + O (e)] 



(A.3) 



45 




e^ + O (e3) 



(A.4) 



(A.6) 



(A.5) 
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A. 2 Two-loop sum-integrals 

We need three two-loop sum-integral that are listed below: 

1 



T 

-'-sun 

P2 - (4/d)p2 



Oie) , 



4(47r)4 Ka-kT) 



1 



1 89 vr 
7 36 2"^ 



10 



7£ + 47I - 871 + 0(e) 



'pqP^Q\P + QY 4(47r)4 UttT/ 



4e 



1 19 

- + — + 47B + 0{e) 
e 6 



(A.7) 

(A.8) 
(A.9) 



The setting-sun sum-integral was first calculated by Arnold and Zhai in Refs. The 
remaining two- loop sum-integrals were calculated by Braaten and Petitgirard |28, 29] using 
the techniques developed in ^. 

A. 3 Three-loop sum-integrals 

We need the following three-loop sum- integrals: 

1 



'-ball 



'PQR P^Q^R\P + Q + Rf 



24(47r)2 V47rr 



6e 



1 91 
e 15 



C(-i) 



C(-3) 



(A.IO) 



ball 



PQR 

P^Qm^{P + Q + Rf 



8(47r) 
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47rT 



6e 



1 



17 



Ci-i] 



+ nlE 17 + 157i5 + 12 



6 

C(-i) 



C(-i 



1 

-+ 
e 



(All) 



and 



1 

pP^ 



mp)r ■ 

rp2 



(47r)2e 

4(47r)4 V47rr/ 
' 1 



n(p) 



1 



4 C' -1 
_ + — ^ + 47F 

3 C(-l) 



+ 



+ 



457r2 



46 - 875; - I67I; - 10471 - 247ij log(27r) + 241og^(27r) + — — + 



C'(-l) C"( 

+ 24 ;; J +2- ^ 



C(-i) C(-i 

where the self-energy n(P) is defined as 



- + 167.^ 



+ Ci + 0{e: 



n(P) = rf 



^qQ2{P + Q)2 ' 



(A.12) 



(A.13) 
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and C^^n 



48.7976 and Cj = —38.4672. The massless basketball sum-integral was first 
calculated by Arnold and Zhai in Refs. The sum-integral Eq. ([A.ll| ) was calculated by 

Gynther et al. in Ref. [IS]. The expression for the sum- integral Eq. ( |A.12|) was calculated 
in Ref. pi. 



A. 4 Four-loop sum-integrals 



We also need a single four-loop sum-integral which was calculated in Ref. [15|: 



(47r)2e 



16(47r) 



- + ( 41og ^ + — + 4^^^ ) - + (21og + 7ij)2+ 



AttT 



C(-i) 



AttT 



^ 2^E + 4^^ ) (2 log -;^+1e)+ Ct'^iangle 



512(47r) 



C(-3) / ' "° AttT 

+ i'jE + Qriangle 



-+81og/- ' ' 



+ 0{e) , 



(A.14) 



where C^^,^^^,^ = -25.7055 and C^^^^,^ = 28.9250. 



B. Three-dimensional integrals 

Dimensional regularization can be used to regularize both the ultraviolet divergences and 
infrared divergences in 3-dimensional integrals over momenta. The spatial dimension is 
generalized to d = 3 — 2e dimensions. Integrals are evaluated at a value of d for which they 
converge and then analytically continued to d = 3. We use the integration measure 



3317 • (B.l) 



Air J J (27r) 



B.l One- loop integrals 

The one-loop integral is given by 



2\n 



Jp {P^ + m 

87r r(2)r(n) 



- 22 - 



Specifically, we need 



h 
h 
h 
h 



log(p^ + m^) 
Gvr \2mJ 



V 9 



^47r \2rn) 

SiTm \2mJ 
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327rm'^ \2mJ 
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, 4 vr^ 



e^ + O fe^ 



(B.3) 
(B.4) 
(B.5) 
(B.6) 
(B.7) 



B.2 Two-loop integrals 



We need the following two-loop integral 
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4(47r)2 \2mJ 
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^ + 6-81og2+ ( 36-y-481og2 + 81og2 2) e + 0{e'^^ 



.(B.8) 



This integral was calculated to order e in Ref. [19| and to order e in Refs. |28, 29| 
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B.3 Three-loop integrals 

We need the following three-loop integrals: 



-^ball 
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The massive basketball was calculated in Ref. |19] to order e'^, and to order e in Ref. 
-^baii obtained by differentiation of /ball with respect to m. The 3-loop integrals J 

and K are calculated in Appendix C. 



B.4 Four-loop integrals 

We need the following two four-loop integrals 
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The triangle diag ram was calculated in Ref. ||30| . The diagram -^[^j-iangic follows from the 
triangle diagram upon differentiation with respect to m?. 
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B.5 Five-loop integrals 
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where = 0.443166. The integrals are calculated in Appendix C. 
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C. Explicit calculations 



In this appendix, we calculate explicitly some of the multi-loop vacuum diagrams in three 
dimensions. 



The three- loop integral J in Eq. ( B.ll ) can be written as 



J = I Kubblc(p)] : 
Ip 



where 



bubble W / (g2+^2)2(p + q): 



^ -|- m-^ 
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By power counting it is easy to see that both J and -^bubble finite in three spatial 
dimension. The latter then reduces to 



-^bubble l^*) 
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Inserting Eq. ( |C.3| ) into Eq. ( C.l ) and using Eq. (|B.5| ) with e = and a mass of 2m. we 
obtain Eq. (|b1ID . 



The integral X can be calculated by noting the relation 

/C,u = -2if-3J. (C.4) 
The integral /rung in ( B.15| ) can be written as 

/rung = / IbnhhleiP) ) (C-5) 

where 

f I 1 

/bubble(/') = / 2? ^ \2~, 2 ■ ^^-^^ 

The integrals /rung and /bubble (p) are convergent in three dimensions. The latter then 
reduces to 

1 p 

huhUeip) = — arctan— . (C.7) 
47rp 2m 

/rung can now be easily found and the result is given by Eq. ( p.l5| ). 
The diagram appearing in J^^i can be written as 

/doublesun = / l~n ; oTo-^sunl/') ) (^'•8) 

Jp ip'^ + m^J^ 

where /sun(p) is 

/■ 1 1 1 

Jqr Q + + (P + q + rj + ?7i 

In order to isolate the divergences in ( |C.8| ), we add and subtract /sun(p = im), and rewrite 
it as 



^doublesun 



1 



X |[4un(p) - /sun(P = + 2Isnn{p)IsuniP = im)- 

-Il^{p = im)] . (C.IO) 

We denote the three terms above by /dsi, /ds2) and /ds3- We first consider /^si. The 
difference /sun (p) — /sun (p = imn) is finite and can be calculated directly in three dimensions. 
We obtain 



/sun(p) - /sun(P = im) = 
1 / /X 



(47r)2 \2m 



3m p 1 , p + 9m , , 

arctan 1 — m — ^ h C(eJ 

p 3m 2 64m^ 



(C.ll) 
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The first term I^si is finite in three dimensions. Using Eq. ( |C.ll ), we obtain 
/dsi = [61og^2 - 61og3 + 41og2 + 5Li,(i) + 



The second term /(js2 can be written as 



Id: 



s2 



2/sun(p = im) 



1 



1 



1 



1 



ipqr {P^ + ^Ti^)^ + + (p + q + r)^ + 

Using Eqs. ( |B.^ ) and ( B.1C| ), we obtain 



(C.12) 



(C.13) 



1 



lOe 



.s2 



16m(47r)5 \2mJ 

4^ + (8 - 12 log 2)- + 52 + ^vr^ - 96 log 2 + 44 log^ 2 + 0{e] 



(C.14) 



Similarly, /dss can be written as 

^ds3 = -4un(P = 
1 



32m(47r)5 V2m/ 



lOe 



1 



1 



TT 



^ + (12 - 161og2)- + 108- — - 1921og2 + 801og^2 + C'(e 



- e 12 
Adding Eqs. (|aT2D , (|Cl§), and (|Cl5| ), we obtain 

1 /j^yoe 



(C.15) 



-^doublesun 



32m(47r)5 \2mJ 

X 



1 131 

^ + (4-81og2)- -4 + — vr^ -961og3 + 641og2 + 1041og^2 + 



12 



+ 80Li2(|) + O(e) 



(C.16) 



Let us finally discuss the five-loop integral appearing in Eq. ( B.17 ). It can be written as 



'4j 



bubble (p) i^triip, q)Y 



pq 



where 



ntri(p, q) 



J.2 _|_ ^2 _|_ + m? (q + r)2 + 



The diagram ( C.18| ) is finite in three dimensions and can be written as [31, 32 1 

arctan(\/D/C) 



ntri(p,g) 



where 



C 



D 



+ + p • q + 4m^ 



p q {p - q) + Am-^lp-^q^ - (p • q)''] 
4^6 



(C.17) 
(C.18) 

(C.19) 

(C.20) 
(C.21) 
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The integral (IoTtI) can now be evaluated numerically by first averaging over angles and 
then integrating over p and q. This yields 



m(47r)5 \2mJ 



lOe 



[0.443166] . 



(C.22) 
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